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Abstract 


The  propagation  of  the  electromagnetic  signals  generated  by 
localized  electromagnetic  current  distributions  in  perfect  crystals  is 
analyzed  under  the  hypotheses  that  dispersion  and  magnetic  anisotropy 
are  negligible  and  all  wavelengths  are  large  compared  with  interatomic 
distances.  The  signal  fields  are  shown  to  converge  asymptotically  in 
energy,  for  t -►  <»,  to  the  sum  of  two  waves  which  propagate  outward  from 
the  source  region  with  the  group  velocities  of  the  crystal.  The 
polarizations  and  waveforms  of  these  waves  are  calculated  explicitly 
in  terms  of  the  crystal  parameters  and  the  source  field. 
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§1.  Introduction. 


This  paper  presents  an  analysis  of  the  propagation  of  electromag- 
netic signals  in  perfect  crystals  under  the  following  physical  hypotheses: 

• The  sources  of  each  signal  are  electric  currents  that  are 
localized  In  a bounded  portion  of  the  crystal. 

• The  frequencies  contained  in  the  sources  are  limited  to  a band  in 
which  dispersion  and  magnetic  anisotropy  in  the  crystal  are 
negligible. 

• The  wavelengths  corresponding  to  the  frequency  band  of  the  sources 
are  large  compared  with  the  interatomic  distances  in  the  crystal. 

Each  such  signal  can  be  characterized  by  its  electric  field  E and  magnetic 
field  H which  satisfy  Maxwell's  equations  with  a tensor  dielectric 
constant  £ and  a scalar  magnetic  permeability  Pg  that  are  characteristic 
of  the  crystal.  Thus 

(1.1)  V X H - £ 8E/9t  = J , V X E + Ug  3H/3t  = 5 

where  J is  the  density  of  the  electric  currents  that  generate  the  signal. 
It  is  assumed  that  J = 0 for  all  t < t^  and  the  signal  generated  by  J is 
characterized  as  the  solution  of  (1.1)  that  satisfies  the  initial 
condition 

(1.2)  E = 0,  H = ^ for  all  t < tg 

A solution  of  the  initial  value  problem  (1.1),  (1.2)  was  given  by 
G.  Herglotz  [3].  An  exposition  of  his  method  may  be  found  in  [2].  This 
solution  makes  it  possible  to  calculate  the  signal  fields  E and  H at  any 


point  in  space-time.  However,  it  does  not  reveal  the  structure  of  the 
signal  or  its  dependence  on  The  analysis  presented  below  is  based  on 
the  author's  work  on  strongly  propagative  anisotropic  media  [8],  The 
main  result  of  that  work  is  a construction  which  assigns  to  each  wave 
field  with  finite  energy  an  asymptotic  wave  field  which  describes  a wave 
that  propagates  outward  with  the  group  speeds  of  the  medium.  The 
difference  between  the  exact  and  as)rmptotic  wave  fields  tends  to  zero  in 
energy  when  t ->•  «>.  The  primary  features  of  the  asymptotic  wave  are 
determined  by  the  medium.  Superposed  on  these  features  is  a waveform 
that  is  characteristic  of  the  wave  sources.  In  this  paper  the  results  of 
[8]  are  used  to  construct  asymptotic  wave  fields  for  electromagnetic 
signals  in  crystals. 

The  remainder  of  the  paper  is  organized  as  follows.  §2  contains  | 

i 

a summary  of  the  properties  of  the  slowness  surfaces  and  wave  surfaces  i 

for  the  Maxwell  system  that  are  needed  in  the  subsequent  analysis.  In  §3  j 

I 

a simole  Hilbert  space  method  is  used  to  solve  the  initial  value  problem  | 

I 

and  derive  the  slowness  surface  representation  of  the  signal  field  E that 
is  the  starting  point  for  the  asymptotic  analysis  of  [8].  In  §A  the 
asymptotic  wave  fields  for  ^ are  constructed  and  their  dependence  of 
their  polarization  and  waveform  on  the  sources  and  crystal  structure  are 
determined. 
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§2.  The  Slowness  and  Wave  Surfaces. 

The  properties  of  the  slowness  and  wave  surfaces  of  general 
strongly  propagative  media  were  discussed  in  [8].  Thorough  discussions 
of  the  case  of  electromagnetic  waves  in  crystals  may  be  found  in  the  books 
of  A.  Sommerfeld  [6]  and  Courant-Hilbert  [2]  where  the  surfaces  are  called 
the  normal  and  ray  surfaces  respectively.  This  section  contains  a review 
of  the  definitions  and  properties  of  these  surfaces  that  are  needed  for 
the  analysis  of  electromagnetic  signal  propagation. 

If  the  field  ^ is  eliminated  between  the  two  equations  (1.1)  the 

equation 

(2.1)  ePo0  + VxVx|  = -po|^ 


is  obtained.  To  solve  this  equation  it  is  convenient  to  use  a system  of 
Cartesian  coordinates  x = (xj,X2>X3)  € with  axes  along  the  principal 
axes  of  the  dielectric  tensor  e.  In  such  systems  the  components  of  e 


have  the  diagonal  form  e 


jk 


The  notation  c^  = (e^  Pj) 
will  be  assumed  that 


-1/2 


where  ej,  £2  and  £3  are  positive, 
will  also  be  used.  For  definiteness  it 


(2.2)  ® ^ ^ 5 *^3 

Equation  (2.1)  is  equivalent  to  the  system 


(2.3) 


at 


+ 

2 J 


a [ail 

axj  [ax. 


, aE,  ^ 3Ea 
ax^  axj 


aj 


at 


, j = 1.2,3 


where  A » a^/3Xj  + a^/ax^  + a^/3Xj  and  E^  and  Jj  are  the  components  of 
and  ? in  the  coordinate  system  (Xj.Xj.Xj).  The  magnetic  fields  corres- 
ponolng  to  solutions  of  (2.3)  can  be  obtained  from  the  second  equation  of 


(1.1)  by  a t-lntegratlon. 
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Plane  waves  in  the  crystal  are  solutions  of  (2.3)  with  = 0 
which  have  the  form 


(2.4) 


Ej  - , j . 1.2,3 


where  the  a^  are  independent  of  t and  x.  Substitution  of  (2,4)  into 
(2.3)  with  = 0 leads  to  three  linear  equations  for  Sj,  a^,  a^  which 
may  be  obtained  by  the  substitutions  9/9t  -►  -io),  iPj  » ^ 

and  Jj  -►  0.  Non-trivial  solutions  exist  if  and  only  if  the  determinant 
of  this  system  is  zero.  This  condition  reduces  to  the  following  disper- 


sion relation  between  the  frequency  o)  and  the  components  p^  of  the  wave 
vector. 


(2.5)  Q(o)^,p)  = o)^(ii)'‘  - Q2(p)iD^  + Q^(p))  = 0 

where 

(2.6)  Qa(p)  = (c|  + c|)  p?  + (c^  + c\)  p\  + (c^  + c|)  p] 

(2.7)  Q4(p)  = |p1^  c|  p2  + c\  c2  p2  + cj  c2  pp 

The  equation  Q(\i,p)  = 0 has  the  roots  0 E P(,(p)  < Ui(p)  < Uz^P)  where 

(2.8)  pj(p)  = j {Q2(p)  - v^rpT}  , yz(p)  = \ {Q2<p)  + . 

(2.9)  D(p)  = Q2(p)^-4QJp)  = [(CiPi+C3P3)^+C^p^][(CjPj  -C3P3)"+C^p^] 
and  Cj,  C2,  C3  > 0 are  defined  by 

(2.10)  - c^3  - c^  , cl  ~ c]  - c\  , C^3  = c"  - c^ 


Hence  the  dispersion  relation  has  the  non-zero  solutions 


r w 


w = A^(p)  = /y^(p)  , a = 1,2 


The  functions  are  homogeneous  of  degree  1.  Hence  to  each  p 0 

there  correspond  two  plane  waves  (2.4)  with  phase  speeds 


(2.12) 


A(j(p)/|p|  = A^(p/|p|)  . a = 1,2 


An  electromagnetic  classification  of  crystals  may  be  based  on  the 
p-dependence  of  the  phase  speeds  (2.12);  cf.  [6].  There  are  three  classes 
which  are  defined  by  the  conditions  (1)  ~ ^3»  (2)  > ^2  “ ^3 

or  Gj  = £2  > £3  and  (3)  £j  > £j  > £3’  Crystals  in  the  first  class  are 
called  isotropic  because  Aj(p)  = Xj(p)  = c|p|  and  the  phase  speeds  are 
Independent  of  p/|p|.  Crystals  in  the  second  class  are  called  uniaxial 
because  Xi(p)  and  X2 coincide  along  a single  axis.  Finally,  crystals 
in  the  third  class  are  called  biaxial  because  Xj(p)  and  Aj (p)  coincide 
along  the  two  distinct  axes  through  the  origin  and  the  points  (/cj-Cj , 0, 
±/c|-c|)  [6].  Note  that  by  (2.9)  D(p)  is  the  square  of  a polynomial  if 
and  only  if  at  least  one  C^  is  zero.  Hence  the  biaxial  crystals  are 
precisely  those  for  which  Q(y,p)  is  irreducible.  In  the  remainder  of 
the  paper  only  the  biaxial  case  is  discussed.  Results  for  the  other 
cases  can  be  obtained  by  allowing  one  or  both  of  the  constants  £j  to 
coincide. 

The  variation  of  the  phase  speeds  (2.12)  with  p/jp]  can  be 
visualized  by  means  of  the  slowness  surface  S = Sj  U defined  by 

(2.13)  * {p  6 R^;  A^(p)  = 1}  = {p  € R^:  |p|  = A^*(p/|p|)},  a = 1,2 


The  definition  of  the  roots  A^(p)  implies  that  S Is  the  real  algebraic 
variety 


J 


6 


r 


(2.14) 


S = {p  e R^:  Q(l,p)  = 0} 


It  is  a member  of  the  class  of  slowness  surfaces  whose  geometry  was 
studied  in  [8],  Each  is  a closed  surface  which  is  star-shaped  with 
respect  to  the  origin.  Although  Q(l,p)  is  irreducible  its  restrictions 
to  the  three  coordinate  planes  split  into  quadratic  factors  and  the 
sections  of  S in  each  of  these  planes  consist  of  a circle  and  an  ellipse. 
It  is  clear  from  (2.8),  (2.9)  that  the  set  Zg  = Sj  n Sj  of  singular 
points  of  S consists  of  the  four  points  in  which  the  two  axes  where 
^i(p)  “ intersect  S.  The  two  sets  - Zg  are  disjoint  analytic 

surfaces.  It  is  known  that  S is  a Fresnel  surface  with  equation  [6] 


(2.15) 


S 


'p  e R 


3 . 


Fresnel  surfaces  were  studied  extensively  in  the  nineteenth  and  early 
twentieth  centuries;  see,  e.g.  [1,  5]. 

It  will  be  shown  in  §4  that  electromagnetic  signals  in  crystals 
propagate  outward  from  the  source  region  with  the  group  velocities 
defined  by 


(2.16) 

The  homogeneity  of  ^^^(p)  and  (2.13)  imply  that  (2.16)  defines  a point  x 

on  the  tangent  plane  x • p = 1 to  at  p.  The  variation  of  the  group 

speeds  v = |V  X I with  the  propagation  direction  0 = x/|x|  = V X /|V  X 1 
f a'pa‘  r f ej  I'pa'pol' 

can  be  visualized  by  means  of  the  wave  surface 

(2.17)  W=»{xeR^:  x*p  = lisa  tangent  plane  to  S} 

W is  the  polar  reciprocal  of  S [2,  8].  It  is  known  that  it  is  also  a 
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Fresnel  surface,  namely  the  surface 


(2.18) 


W = /x  e R^:  I 


= 0 


j=l  c^^x|^-l 


obtained  from  (2.15)  by  the  substitutions  p ->•  x,  c^  -►  c^  \ In  particular, 
W = Wj  U Wj  where 


(2.19)  W = {x  G R^:  X’(x)  = 1}  , a = 1,2 

m ct 

j 

and  X' (x)  is  obtained  from  X (p)  by  the  same  substitutions.  It  follows 

a a ^ ^ 

that  there  are  exactly  two  group  velocities  in  each  direction  6 = x/|xl 
except  the  directions  of  the  four  singular  points  of  W.  The  corresponding 
group  speeds  are  given  by 

(2.20)  v^(9)  = 1/X^(0)  , a = 1,2 

It  is  clear  from  (2.8)  that  0 < V2(0)  < v^(0)  for  all  0. 

The  Gauss  map  N:  S = {0  e R^:  |0|  = 1}  and  the  polar 

reciprocal  map  T:  S W are  also  needed  in  §4.  By  definition,  N assigns 

to  each  p G S the  set  of  all  exterior  unit  normal  vectors  to  S at  p;  see 

[8]  for  the  construction  of  N(p)  at  singular  points.  In  the  case  of  the  i 

Fresnel  surface  (2.15)  N is  single-valued  and  analytic  except  at  the  set 
Zg  of  four  singular  points,  T is  defined  by 

(2.21)  T(p)  = (p.N(p))"^  N(p)  , pG  S 


It  is  known  that  p • N(p)  is  bounded  away  from  zero  [8].  Hence  T is  also 
single-valued  and  analytic  on  S - Zg.  To  each  p G Zg  it  assigns  one  of 
the  four  circles  defined  by  the  intersections  of  W with  the  four  planes 
/c^-c j Xj  ± /c|-c|  Xj  ± Cj/c^-c^  = 0.  The  union  of  these  circles  will 


I 

t 
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be  denoted  by  = TZg.  Note  that  x = T(p)  is  given  by  (2.16)  at  the 
regular  points  of  S.  Since  S is  also  the  polar  reciprocal  of  W the  map 
T S W -*■  S has  the  same  form  as  T;  that  is,  T ^ (x)  = (x  • N'(x))  * N'(x) 
where  N':  W is  the  Gauss  map  of  W.  T ^ is  single-valued  and 

analytic  on  W - Z^  and  T ^Z^  = Zg  is  the  set  of  four  circles  in  which 
S intersects  the  four  planes  /c“^-c“^  Pj  ± /c“^-c“^  p^ + c^  ^ /c'^^-Cj^  = 0. 
It  was  shown  in  [8]  that  if  Zg  = Zg  U Zg  and  Z^  = Z^  U z^  then  T is 
bijective  and  analytic  between  S - Zg  and  W - Z^, 

Finally,  the  inverse  Gauss  map  N -»■  S is  needed  in  §4.  It 

is  clear  from  (2.21)  and  the  geometry  of  W that  N ' is  double- valued  for 
all  6 € - Zg  where  Zg  is  the  intersection  of  with  the  four  rays 

through  the  singular  points  of  W.  More  precisely,  if 

(2.22)  x^“^(9)  = v^(0)  0 , a = 1,2 

then 


(2.23) 


p(“)(0)  = x^“^(0)  , a = 1,2 


defines  the  two  branches  of  N * and  each  p^*^^  is  analytic  on  - Zg 
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§3.  Representation  of  the  Signals  by  Slowness  Surface  Integrals. 

In  this  section  Fourier  analysis  Is  used  to  represent  the  signals 
as  slovmess  surface  Integrals.  The  analysis  follows  the  method  of  [8,  9] 
to  which  reference  Is  made  for  analytical  details.  Only  the  main  steps 
are  given  here. 

The  Maxwell  equations  (2.3)  may  be  written  as  a matrix  system 


(3.1) 


3"u 

¥t^ 


+ A u = f (t,x) 


where  u = u(t,x)  = (Ej,E2»E3)  (T  denotes  the  transpose),  f = -e  9J/3t, 

T 

e = dlag  (ej,£2,e3)  and  J = (Jj, 32,33)  . A Is  the  3x3  matrix 

differential  operator  defined  by  A = c^  A^ , c^  = dlag  (Cj.c^.cp  and 

Aq  = (9^/9x.9x,  - 6., A).  The  solution  of  the  Initial  value  problem  for 
3 R 3 R 

(3.1)  will  be  constructed  In  the  Hilbert  space  3f  = L2(R^,C^)  with  scalar 
product 


(3.2)  (u,v).^  = f u(x)*  c ^ v(x)  dx 

Jr3 

where  u = u Is  the  Hermltlan  adjoint  of  u.  If  u G JC  then  Au  Is  well- 
defined  as  a distribution  and  a linear  operator  A:  Uf  -»■  If  Is  defined  by 


(3.3) 


D (A)  = 3fn{u:  AuGUC},  Au  = Au 


It  can  be  shown  by  Fourier  analysis  that  A Is  selfadjolnt  and  non-negative: 
A*  = A > 0. 

System  (3.1)  will  be  Interpreted  as  an  ordinary  differential 

equation 

(3.4)  ^ + Au  = f(t,.),  t e R 

for  a function  t -*•  u(t,*)  G Tf.  To  Integrate  (3.4)  let  be  the 
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1/2 

non-negative  square  root  of  A and  note  that  D(A  ) is  a Hilbert  space 


, 1/2 


under  the  graph  norm  of  A . If  it  is  assumed  that  the  source  function 


J € (R,D(A^^^))  and  supp  J C {t:  Iq  < t < t^  + X}  then  f = -e  dJ/dt 

e C(R,D(A^/^)) , supp  f C {t:  t„  < t < t.  + T}  and 


(3.5) 


1: 


to+T 


dt  = 0 


The  signal  generated  by  f is  by  definition  the  solution  of  (3.4)  that 
satisfies  u(t,*)  = 0 for  t < t^ . It  is  given  by  the  Duhamel  integral 


(3.6)  u(t,*)  = 


{cos  (t-T)  I f(t',*)  dt'  dx,  t > tj 


Indeed,  the  hypotheses  on  f and  the  spectral  theorem  imply  that  (3.6) 
defines  a function  uG  C^(R,3f)  n (R,D(A^^^) ) C(R,D(A))  that  satisfies 

(3.4).  Note  that  since  J and  f have  real-valued  components,  u(t,x) 

= Re  (v(t,x)}  where 


(3.7)  v(t,-)  = e [ f(t',*)  dt'  dx 

•'tn  ■'t„ 


Moreover,  (3.5)  implies  that  for  all  t > t^  + T 


-ItA^/^ 

(3.8)  v(t, •)  = e h,  h 


^ ^1/2  rT 

j e^^^  f(t',-)  dt'  dx 

bn  tft 


The  wave  function  (3.8)  will  be  constructed  by  Fourier  analysis. 
To  this  end  let  <t:  JC  -*■  Jf  denote  the  unitary  operator  defined  by  the 

Fourier  transform.  Then  A = <1>  A(*)'J>  where  A(*)  denotes  the  matrix  multi- 
plication operator  defined  by  A(p)  =c^A(|(p),  A (p)  = (|p|^6..-p.  p.), 

J K .1  ^ 


p e R^.  To  complete  the  spectral  analysis  of  A(p)  note  that  its  char- 
acteristic polynomial  is  the  polynomial  Q(y,p)  with  roots  Uj(p) 

= 0,Pi  (p)  .VijCp)  defined  by  (2.8).  Thus 


(3.9) 


(U  1 - A(p))"'  = I (U  - U„(p))“  P„(P) 

a=0 


where  the  matrices  P^(p)  are  given  by 


(3.10) 


p = cof(yl-A(p)) 

9Q(y,p)/3y 


, a = 0,1,2 


y=yc^(p) 


and  cof  M denotes  the  cofactor  matrix  of  M.  It  follows  as  in  [8]  that 

^a=0 


r(t,‘)  = exp  (-ltA^/^)h  = V^(t,*)  where 


(3.11) 


v^(t,x)  = 3^2 


a 


i(x*p-tX  (p))  ^ 

(2,)».  J , ' "(P'  “P 

K 


X^(p)  = /y^(p)  as  in  §2  and  h = $h.  In  particular,  Vo(t,x)  = v„(x)  is 
the  static  component  of  v(t,x).  The  remaining  two  integrals  may  be 
combined  as  in  [8, §4].  The  result  is  the  representation 

.OO 

(3.12)  v(t,x)  = Vj(x)  + 1 e w(Xx,X)  X^  dX 

^ 0 

where 

(3.13)  w(x,X)  = (2Tr)"®^2  [ e^*'P  P(p)  h(Xp)  |T(p)r*  dS 

Js 

A ^ 

is  an  integral  over  the  slowness  surface  S = Sj  u S,  , P(p)  = Pqj(p)  for 
p g and  T is  the  jwlar  reciprocal  map  of  §2.  Tins  representation  is 
the  starting  point  for  the  analysis  in  [8]  of  the  asymptotic  behavior 
of  v(t, •)  in  K. 


■ • .«  'Mr:-: 
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§4 . Electromagnetic  Signal  Structure  in  Crystals. 

Consider  a source  function  J(t,x)  which  is  localized  in  both 
space  and  time: 

(4.1)  supp  J c {(c,x):  < t < tj  + T and  |x-Xp|  < 6^} 

The  signal  generated  by  J will  reach  points  x in  the  far  field  of  the 
sources,  characterized  by  |x-Xj|  » 6g , after  a time  interval  of 
magnitude  comparable  with  Ci|x-X(j|.  Hence  the  far  field  form  of  u(t,x) 
coincides  with  its  form  for  large  t.  The  latter  is  given  by  the  asymptotic 
wave  functions  of  [8]. 

The  analysis  of  [8],  applied  to  (3.12),  (3.13),  implies  that  the 
signal  u(t,x)  = Re  {v(t,x)}  generated  by  f = -E  ^ 3j/3t  satisfies 

(4.2)  11m  llu(t,*)  - Ug  - u (t,*)llTf  = 0 

where  Ug (x)  = Re  {vp (x)}  is  the  static  component  of  the  signal,  defined 
00 

in  §3,  and  u (t,x)  is  the  asymptotic  wave  function  for  u.  The  latter 
has  the  form 

2 

(4.3)  u“(t,x)  = |x|~’  I G(x  * p^°^^  (0)  - t,p^“^  (0)) , x=  |x|0 

a=l 

(a) 

where  p , a = 1,2,  are  the  branches  of  the  inverse  Gauss  map  of  S 
uefined  in  §2.  Note  that  T p^“^(0)  = v^“^(0)0,  by  (2.22),  (2.23),  and 
p • T(p)  = 1 for  all  p e S.  Hence 

(4.4)  x'p^^^^d)  = lxl0*p^“^(0)  = lx|/v^(0) 

and  (4.3)  represents  the  sum  of  two  waves  which  propagate  outward  with 
the  group  speeds  The  form  (4.3)  of  the  wave  is  determined  solely 


t- 
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by  the  crystal.  The  shape  of  the  wave  profile  G:  R x s R*  depends  on 

both  the  crystal  structure  and  the  source  function  f.  It  is  given  by 
G(T,p)  = Re  {F(t,p)}  where  F:  R x s C®  is  defined  by 


(4.5)  F(t,p)  = n(p)  [(2Tr)"'^"  (|<(p)  | 


e^^^  h(Xp)  X dX 


and 


(4.6) 


n(p)  = 1k(p)1  lT(p)r*  P(p) 


where  K(p)  is  the  Gaussian  curvature  of  S at  p,  \J)(p)  = exp{i^  (v  (p)-v  (p))} 
and  V^(p)  and  V~(p)  denote  the  number  of  principal  ''urvatures  of  S at  p 
that  are  positive  and  negative,  respectively.  Note  that  II(p)  and  i|>(p) 
are  defined  and  analytic  for  all  p € S ~ Zg.  The  integral  in  (4.4) 
converges  in  the  Hilbert  space  3C(S)  with  norm  defined  by 


(4.7) 


*^*!R(S) 


<00  i 

= J J F(x,p)*  c"^  F(t,p)  1k(p)  T(p)|  dS  dx 


and  the  map  0:  3C  ->■  3C(S)  defined  by  0h  = F is  an  isometry  [8] . 

The  profile  (4.5)  corresponding  to  the  signal  (3.8)  can  now  be 
calculated.  Note  that 
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Finally,  integration  by  parts  with  respect  to  t and  use  of  (3.5)  gives 


(4.10)  P(p)  h(Xp)  = (2TT)‘/2i  X"  P(p)  f(-X,Xp),  p e S - Z, 


where  f(co,p)  denotes  the  four-dimensional  Fourier  transform  of  f(t,x). 
Thus,  since  II(p)  has  real  components,  the  profile  G(x,p)  = Re  {F(t,p)1 
generated  by  f has  the  form 


(4.11) 


G(x,p)  = n(p)  Re 


|i>l^(p)  f 


e^”^^  f(-X,Xp)  dX 


Note  that,  since  f(u),p)  is  the  Fourier  transform  of  a real-valued  function 


(4.12) 


f(-X,Xp)  dX  = Sj(-X,p)  - i a^(-x,p) 


where  s^(x,p)  and  o^(x,p)  are  a pair  of  Hilbert  transforms  with  respect 
to  X [7].  Moreover,  an  application  of  Fubini's  theorem  gives 


(4.13)  Sj(x 


, 1 
.P)  = 2 


f(-X,Xp)  dX 


= 1- 
4tt  J, 


f (p  • x+  x,x)  dx 


Hence  the  signal  waveform  s^(x,p)  is  the  integral  of  f(t,x)  over  the 
space-time  hyperplane  t = p*  x+x.  Thus  s^(x,p)  is  essentially  the  four- 
dimensional Radon  transform  of  f(t,x)  [4]. 

Finally,  note  that  the  phase  factor  i/^(p)  in  (4.10)  satisfies 
'i^(p)  = “i  at  elliptic  points  of  S and  (J^(p)  = 1 at  hyperbolic  points. 
Moreover  S^,  the  inner  sheet  of  S,  is  convex  and  has  only  elliptic  points 
while  Sj,  the  outer  sheet,  is  divided  into  five  components  by  the  four 
circles  which  make  up  the  set  Zg  of  parabolic  points.  The  union  of  the 
four  components  that  contain  singular  points  is  the  set  S of  hyperbolic 
points  of  S and  S = s"*"  u S u Z„,  a disjoint  union.  Combining  (4.11), 


I 

t 

I 
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(4.12),  (4.13)  and  these  remarks  gives  the  signal  waveform 


j 

^ n(p)  s^(-T,p), 

p 6 s 

1 (4.14) 

G(x,p)  = i 

1 

( n(p)  a^(-x,p), 

p € s 

where  s^(T,p),  a^(x,p)  is  the  Hilbert  transform  pair  defined  by  (4.12). 


1 
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